Anisotropy and internal field distribution of MgB 2 in the mixed state 

at low temperatures 
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Magnetization and muon spin relaxation on MgB2 were measured as a function of field at 2K. 
Both indicate an inverse-squared penetration depth strongly decreasing with increasing field H 
below about f T. Magnetization also suggests the anisotropy of the penetration depth to increase 
with increasing H , interpolating between a low H c i and a high H C 2 anisotropy. Torque vs angle 
measurements are in agreement with this finding, while also ruling out drastic differences between 
the mixed state anisotropics of the two basic length scales penetration depth and coherence length. 

PACS numbers: 74.25.Op, 74.20.De, 74.25. Ha, 74.70.Ad 



The understanding of the physical properties of the re- 
cently discovered 39 K superconductor^. MgB 2 has made 
rapid progress in the last 3 yearsi 2 - A central issue of 
research has been the involvement in superconductivity 
of two sets of bands with different dimensionality and 
pairing strength£^& This "two-band superconductivity" 
leads to an array of unusual superconducting properties 
such as specific heat^ particularly also to a very unusual 
behavior of the superconducting anisotropics. 6 

For example, a pronounced temperature T dependence 
of the anisotropy jh of the upper critical field H C 2, di- 
rectly related to the coherence length £, was observed^ 
and calculated based on the two-band modelfSa 1 . Strik- 
ingly, calculations of the low field penetration depth 
anisotropy 7a , predicted a much lower anisotropy of this 
quantity, with a T dependence opposite to the one of 
H C 2^ This was experimentally confirmed as well, based 
on measurements, e.g., of iJ c xji2ii& by small angle neu- 
tron scattering (SANS)i^*i^ and scanning tunneling spec- 
troscopy (STS)ii 6 " However, e.g. the experiment of Ref. 
IT5I indicates that whereas in the limit of very low fields 
H 7a is indeed close to 1, it is rising with increasing H, 
as deduced earlier more indirectly 7 

The behavior of the anisotropics of the length scales 
in the mixed state H c \ < H < H C 2 still needs to be clari- 
fied. One point of view^£ surmises constant (with respect 
to H) anisotropies of the penetration depth 7a and the 
coherence length jh, which are, however, different from 
each other. This difference was predicted to lead to a sign 
reversal in the angle dependent torque>i£ Another point 
of view is that these anisotropies are not drastically dif- 
ferent from each other, but both increase with increasing 
field, interpolating from the H c i anisotropy in low fields 
to the H C 2 anisotropy in high fields^ 

Here, we support the latter point of view by analyzing 
SQUID (superconducting quantum interference device) 
and torque magnetization data measured on a MgB2 sin- 



gle crystal with very low pinning, and muon spin relax- 
ation (/iSR) data measured on randomly aligned MgB2 
powder. In the absence of a more elaborate model of the 
mixed state of a two-band superconductor ^» we base the 
analysis on the London model, allowing however for a 
H dependent penetration depth, which is obtained from 
the H or angle 8 dependence of the bulk magnetization 
(SQUID/torque) as well as from the average variation 
of the internal field (/iSR). From SQUID and fiSR we 
find a rapid decrease of the inverse-squared penetration 
depth 1/A 2 ("superfluid density") with fj, H increasing 
below about 1 T, and SQUID and torque data agree on 
the anisotropy 7a increasing strongly with H. The anal- 
ysis of the torque data further suggests that is not 
very different from 7a. 

Single crystals of MgB2 were grown with a high pres- 
sure cubic anvil technique^ and a crystal with particu- 
larly low pinning was selected for measurements with a 
Quantum Design MPMS-XL SQUID magnetometer and 
a noncommercial torque magnetometer^ The crystal 
has a wedge shape, with one of the faces parallel to the 
ab planes. The /xSR experiment on poly-crystalline MgB2 
was similar to the one of Ref. l2ll 

The magnetization of the single crystal was measured 
as a function of H; as can be seen in Fig. ^ the irre- 
versibility is very low above about 0.15 T. Larger irre- 
versibility in lower fields may be due to geometrical bar- 
riers, which is why we did not attempt to directly extract 
H c \. The curves shown are not corrected for demagne- 
tizing effects (the exact demagnetizing factor is difficult 
to estimate due to the sample shape), but we verified 
that any reasonable demagnetization correction does not 
noticeably affect above 0.2 T the results discussed below. 

Within the London model of a standard supercon- 
ductor, the magnetization is proportional to the loga- 
rithm of the applied field, not too close to either H cl or 
H C 2 . Keeping to an analysis within the London approach, 
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FIG. 1: Magnetization M vs field fj, a H at 2 K on a MgB 2 
single crystal with H||c and H\\ab. 



but dropping the requirement of a constant penetration 
depth A, we have 1/A 2 oc dM rev /dlni?. The so obtained 
1/A 2 is plotted in Fig. To avoid overloading the graph, 
only curves assuming M rcv = (M#f + Mh\,)/2 are shown; 
except in the lowest H using instead Mh] or Mhi leads 
to very similar results. The curves for both H directions 
were normalized by the same constant factor. The shaded 
box indicates the low H region, where we are uncertain 
about the obtained penetration depth because of i) irre- 
versibility, ii) demagnetizing effects, and iii) deviations 
from the London model due to the vicinity of H c \ (see 
below). 

Also plotted in Fig. [2 is 1/A 2 obtained from the muon 
spin depolarization rate a measured on randomly aligned 
powder at the same temperature. The depolarization 
rate a is a measure of the average variation of the in- 
ternal field within a superconductor, and in the mixed 
state (again not too close to one of the critical fields) is 
directly proportional to 1/A 2 , since A is the fundamental 
length scale of the variation of the field in a superconduc- 
tor (cf. Ref . l2lh . An issue to be aware of when deducing 
A in this way is the possible influence of pinning, which 
can lead to an extrinsic increase in er. In a previous /iSR 
experiment on MgB2, the whole H dependence of a was 
indeed ascribed to pinningi 2 ^ To check for the possible 
influence of pinning on a (as opposed to the magnetiza- 
tion), we performed time-dependent measurements of a 
in several fields: After reaching 2K (field cooled), statis- 
tics was gathered for 10 min, then stopped and restarted 
(repeated 5 to 10 times). Except for 0.1 T, changes of a 
with time are well below error bars, and no clear trend 
discernible. This suggests that for higher H even at 2K, 
pinning is not influencing cr much, and the H dependence 
of a indeed intrinsic. That we observed a very similar H 
dependence of a in samples from two sources synthesized 
slightly different supports this; an intrinsic cr(H) depen- 
dence was also proposed in Ref. f23t 

Concerning the dM/d(lnff) curves, it may be argued 
that the low field behavior is not unexpected even for 
a normal superconductor, since in the limit H — > H c \ 
it is expected 2 ! dM/d(lnff) oc H/(H-H cl ). However, 
dM/d(ln-ff) should reach values close to the normal Lon- 
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FIG. 2: Comparison of 1 /A 2 vs H obtained from dM/d In H of 
Fig. Q and from the muon spin depolarization rate measured 
on unaligned powder (circles). The shaded box indicates fields 
close to or lower than H c i (see text). 



don slope rather quickly (within 2 — 3H c i) and the vari- 
ation presented in Fig. |3 is spread over a considerably 
larger field range. Furthermore, the influence of the 
vicinity of H c i on a is opposite. The close similarity 
of the H dependence of the penetration depth obtained 
from rather different quantities (bulk magnetization from 
SQUID and internal field variation from /iSR) strongly 
suggests that all curves in Fig. [5] indeed show the X(H) 
dependence, outside of the shaded box indicating vicinity 
of H c \. A similar strong depression of 1/A 2 with H was 
deduced previously (for H\\c) from an analysis of SANS 
form factors^ The SQUID curves additionally indicate 
that above lOkOe 1/A 2 no longer varies strongly. 

The SQUID measurements in the two field configura- 
tions also yield the anisotropy of A. For H\\c the screen- 
ing currents flow within the ab plane, giving 1/A 2 fc . For 
H \\ab the currents flow also perpendicular to the planes, 
giving l/(A a (,A c ). The ratio of the dM/d In H curves for 
the two field configurations thus corresponds to y\. Con- 
sidering the curves in Fig. |21 we can see that i) in low H 
the anisotropy is very small, ii) in high H j\ is of the 
order of about 6 or even 7, and iii) the variation with 
H of 7a is most pronounced in low H. We stress the 
fact that when considering the high field region alone 
the standard London model with constant A and £ and 
a constant common anisotropy 7 describes the data rea- 
sonably well. This indicates that at low T in high H 
MgB 2 is close to a "standard superconductor" with high 
anisotropy. Corresponding to this is the absence of an 
unusual H C 2(9) dependence at low T, in contrast to the 
situation closer to J , (j i2i2^ 2 & [of course the 7a analysis 

breaks down as fi a H^ ^ H^ (— 2.8 T for this crystal)]. 

An alternative method to determine the penetration 
depth anisotropy is to analyze the angular 9 dependence 
of the torque r in fixed H. We previously used this 
method at much higher T, finding also indications of an 
anisotropy increasing with Hi However, thermal flue- 
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FIG. 3: Angle dependence of torque r in 0.2 T at 8K 
(symbols). Dashed line: theoretical description— assuming 
7a -C 7h; full line: description with 7a — ^h- 



tuations and additional intermixture of the two sets of 
bands by thermally excited quasiparticles, complicate the 
analysis there. To provide a direct comparison with the 
SQUID results and give a quantitative estimate of 7a we 
measured t(0) at low temperature. 

Due to increased irreversibility at lower T, it is impor- 
tant not only to use a crystal with low pinning, but also 
employ the "shaking technique" developed by Willcmin 
et alOL Since the magnetometer equipped with this tech- 
nique cannot reach 2 K measurements were conducted 
at 8, II and 15 K. These temperatures should be low 
enough to avoid too strong an influence of thermal fluc- 
tuations/excitations, as well as to probe the low temper- 
ature limit of the calculated**^ anisotropies. 

The data were analyzed with Eq. (18) of Ref.lT^. allow- 
ing for a difference in the anisotropies of the penetration 
depth and the coherence length. Such a difference was, 
however, not found in any of the curves analyzed, and a 
sign reversal of the torque, a key prediction of Ref. 17 for 
7a <S IHi was never observed [for an example see Fig. 
p3- A preliminary report on this issue is given in Ref. 
HB. The best descriptions with Eq. (18) of Ref . E3 of the 
data were rather achieved for 7a ~ jh- The large num- 
ber of parameters involved and the numerical condition 
of the fit formula result in extended error bars though, 
so that small differences between 7a and cannot be 
completely ruled out (but large differences*! as calculated 
with H independent anisotropies can). 

The resulting field dependence of the anisotropy 7 = 
7a ~ 7h is shown in Fig. 0] The anisotropy is monoton- 
ically increasing with increasing field, up to 1.5 T (the 
maximum attainable by the magnetometer used). In the 
lowest fields, this rise seems to be much steeper than 
above 0.5 T. However, as a cautionary note, even with 
"shaking", the irreversibility cannot be said to be negli- 
gible below 0.1 T and it should also be kept in mind that 
we are approaching H cl (see Fig. 0). 

A report by another group of torque vs angle mea- 
surements performed at 10 K claimed a field-independent 
anisotropy of the order of 7 ~ 4.3r 28 However, analyzing 
the same data, a different conclusion of an anisotropy 
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FIG. 4: Anisotropy determined from an analysis of the torque 
data (Fig. with Eq. (18) of Ref.ll7l at various low temper- 
atures, as a function of field H . Also shown is the anisotropy 
determined from Fig. [5] (o) and the anisotropy very close to 
T c (dotted line, see text). 



that does increase with H, in not too large fields, may 
also be reached m& The results of Ref. |2^, as well as Ref. 
l3fl (not finding a field-dependence as well) may be recon- 
ciled with the ones of Ref. and the present results by 
assuming a tendency of 7 to saturate in higher fields. 

For comparison, a rough estimate of "f(H) from SQUID 
magnetometry in fields along the principal axes [Fig. 
is plotted in Fig. 21 as well. Qualitatively ('y(H) being 
an increasing function) this is consistent with the torque 
results. We attribute the numerical discrepancy to the 
large scattering as visible in Fig.|3and the corresponding 
uncertainty in the estimation of 7. The low T behavior 
of j(H) is in strong contrast to the one very close to 
T c , where between H c \ and H C 2 7 ~ 2 is constant^ as 
indicated by the dotted line in Fig. 2] 

A field dependent anisotropy at low T (2 K) had been 
deduced based on different experiments as well. Bou- 
quet et alml reported a H dependent effective anisotropy 
based on specific heat measurements. Since these are 
sensitive mainly to the coherence length, the experiment 
suggests the anisotropy to be H dependent. Cubitt 
et a/p*£ observed the anisotropy of the vortex lattice jvl 
to increase strongly with increasing H, from less than 
1.5 in 0.2 T to about 3.8 in 0.5 T, Jvl(H) being more 
steep in higher H . Keeping to the London model, the 
anisotropy of the vortex lattice should be equal to the 
penetration depth anisotropy 7a. - 3 Our results extend 
to higher H and agree qualitatively with those of Cu- 
bitt et at, but we do not find a particularly steep ^(H) 
around 0.5 T, but rather a slower field dependence. Very 
recently, Lyard et alM proposed a similar H dependence 
of anisotropies, based on a London analysis of magneti- 
zation data measured at much higher T. 

The strong field dependence at low T of all anisotropies 
obtained from the measurements presented here, as well 
as by other groups, are readily explained qualitatively by 
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a faster suppression with H of superconductivity in the 
more isotropic n bands, increasing the relative contribu- 
tion of the highly anisotropic a bands. Such a suppres- 
sion, consistent with the overall decrease of 1/A 2 (Fig. 
[3J|, was also observed, e.g., by spectroscopic means 
and is not unexpected due to the smaller gap in the tt 
bands. Corresponding larger vortex cores and a lower 
a H^ 2 " have been conjectured from STS and specific heat 
measurements^^ If n and a bands were independent, 
a real upper critical field fM H^ 2 ~ 0.5 T would mark the 
destruction of superconductivity in the 7r bands due to 
vortex core overlap. Since, however, the bands are cou- 
pled together even at zero T, U H£ 2 " degenerates into a 
broad crossover (completely blurred for T — » T c ). Our 
results indicate that this crossover region is very broad, 
extending down to almost zero field. In high H super- 
conductivity in the 7r bands is still induced from the a 
bands likely up to the bulk H c2 , but with a much de- 
pressed order parameter. It should be noted that within 
this picture, there would in principle be two different co- 
herence lengths to consider^ and that in the H region of 
interest, the vortex cores in the tt band overlap^ enough 
to seriously question the applicability of a London anal- 
ysis. This may explain the remaining discrepancies be- 



tween the anisotropies obtained from different measure- 
ments and calls for further theoretical work, although in 
a qualitative way the London analysis works out remark- 
ably well, particularly in high H at low T, in terms of a 
"standard anisotropic" a band only superconductor. 

In conclusion, /iSR and magnetization data show the 
"superfluid density" 1/A 2 in MgB 2 at 2 K to strongly de- 
crease with increasing field below about 1 T. In parallel, 
the penetration depth anisotropy increases, and is not 
drastically smaller than the coherence length anisotropy 
(in the same field). This behavior is due to the fast sup- 
pression of the contribution to superconductivity of the 
more isotropic n bands with weaker superconductivity. 
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